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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. Massive stars, supernovae (SNe), and long-duration gamma-ray bursts (GRBs) have 
a huge impact on their environment. Despite their importance, a comprehensive knowledge of which massive 
stars produce which SN/GRB is hitherto lacking. We present a brief overview about our knowledge of mass 
loss in the Hcrtzsprung-Russell Diagram (HRD) covering evolutionary phases of the OB main sequence, the 
unstable Luminous Blue Variable (LBV) stage, and the Wolf-Rayet (WR) phase. Despite the fact that metals 
produced by "self-enrichment" in WR atmospheres exceed the initial host galaxy - mctallicity, by orders 
of magnitude, a particularly strong dependence of the mass-loss rate on the initial metallicity is found for 
WR stars at sub-solar metallicities (1/10 - 1/100 solar). This provides a significant boost to the collapsar 
model for GRBs, as it may present a viable mechanism to prevent the loss of angular momentum by stellar 
winds at low metallicity, whilst strong Galactic WR winds may inhibit GRBs occurring at solar metallicities. 
Furthermore, we discuss recently reported quasi-sinusoidal modulations in the radio lightcurves of SNe 2001ig 
and 2003bg. We show that both the sinusoidal behaviour and the recurrence timescale of these modulations 
are consistent with the predicted mass-loss behaviour of LBVs, and we suggest LBVs may be the progenitors 
of some core-collapse SNe. 

ABSTRACT 

Massive stars, supernovae (SNe), and long-duration gamma-ray bursts (GRBs) have a huge impact on their 
environment. Despite their importance, a comprehensive knowledge of which massive stars produce which 
SN/GRB is hitherto lacking. We present a brief overview about our knowledge of mass loss in the Hertzsprung- 
Russell Diagram (HRD) covering evolutionary phases of the OB main sequence, the unstable Luminous Blue 
Variable (LBV) stage, and the Wolf-Rayet (WR) phase. Despite the fact that metals produced by "self- 
enrichment" in WR atmospheres exceed the initial host galaxy metallicity, by orders of magnitude, a 
particularly strong dependence of the mass-loss rate on the initial metallicity is found for WR stars at sub- 
solar metallicities (1/10 - 1/100 solar). This provides a significant boost to the collapsar model for GRBs, as it 
may present a viable mechanism to prevent the loss of angular momentum by stellar winds at low metallicity, 
whilst strong Galactic WR winds may inhibit GRBs occurring at solar metallicities. Furthermore, we discuss 
recently reported quasi-sinusoidal modulations in the radio lightcurves of SNe 2001ig and 2003bg. We show 
that both the sinusoidal behaviour and the recurrence timescale of these modulations are consistent with the 
predicted mass-loss behaviour of LBVs, and we suggest LBVs may be the progenitors of some core-collapse 
SNe. 
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1. INTRODUCTION 

Massive star winds and core-collapse supernovae 
(SNe) have a huge influence on their environments 
by driving the chemical evolution of galaxies and 
shaping the interstellar medium over all cosmolog- 
ical epochs, since the very first stars came into exis- 



tence. Despite their importance, the lives and deaths 
of massive stars are p oorly understood. De spite the- 



oretical progress (e.g. Hirschi et al. 2004h . it is not 
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known with any degree of certainty which massive 
stars produce which SNe/GRB. 

While progress is being made in the direct 
identification of SN progenitors by searching for 



these in pre-explosio n images (e.g. ISmartt 2002t 
van Dvk et al. Il2003h . current progenitor masses ap- 
pear to be limited to stellar masses not significantly 
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greater than ~10-15Mq, likely as a result of the ini- 
tial mass function. 

The evolution of more massive stars (M > 40 
Mq) is largely unconstrained, but it is generally ac- 
cepted that mass loss drives these objects through 
the O star, Luminous Blue Var iable (LBV), and 
Wolf -Rayet (WR) phases (e.g. IChiosi fc Maeder 
19861 ). Mass loss also determines the stellar mass 



before collapse, and is therefore relevant for the type 
of compact remnant that is left behind (i.e. neutron 
star or black hole). This process is expected to de- 
pend on_th^_jnetaj_cont«it (Z) of the host galaxy 
(e.g. lEldridee fc Vink I [2006h . As WR sta rs are the 
likely progenitors of long-duration GRBs ( Wooslev"! 
1993 ) , the strength of WR winds as a function of Z 
is especially relevant for setting the threshold Z for 
forming GRBs. 

Furthermore, massive stars explode in environ- 
ments that have been modified by mass loss from 
the progenitor. The SN ejecta interact first with this 
circumstcllar material before interacting with inter- 
stellar material. We might therefore expect that the 
different wind properties over the lifespan of a mas- 
sive star be imprinted onto the resulting circumstcl- 
lar media (CSM), and we would expect these dif- 
ferences to be seen in the interaction between the 
SN ejecta and surrounding material. By quantifying 
these differences one may be able to constrain the 
evolutionary phase of the exploding object. 

Over the last decades, radio observations of SNe 
have provided a means with which to constrain 
the density of the CSM around core-collapse SNe. 
The inferred mass-loss rates from modelling of most 
radio SN light curves yield values of M ^1CP 6 
4 M Q yr 



(e.g. compilation in IWeiler et al 



20021) . Unfortunately, these average mass-loss rates 
are generally only accurate to within a factor of 
~10, and are typical of almost all types of massive 
star, making it difficult to pin down the evolutionary 
phase during which core-collapse occurred. 

However, a small subset of radio SNe show quasi- 
periodic modulations in their radio lightcurves. We 
argue that this type of modulation may be the re- 
sult of an LBV that underwent S Doradus varia- 
tions which entailed opacity changes in the wind- 
driving region, resulting in varying mass-loss rates 



(IVink fc de Koter 1120021: Ikotak fc Vink 1120061) 



Given the crucial role that mass loss plays for 
massive star evolution, we briefly discuss recent 
mass-loss predictions in order of decreasing tem- 
perature: WR stars — > OB supergiants — ► LBVs 
(Sects. 2-4). In Sects. 5 and 6, we link our knowl- 
edge of mass loss to certain types of radio SNe, and 
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Fig. 1. Mass loss versus initial Z for late- type nitrogen- 
rich (WN) stars (solid line) and carbon-rich WC stars 
(dashed line). Note that metal self-enrichment is ac- 
cou nted for, but does not en ter in our expression of Z. 
See I Vink fc de Koter I (|2005l ) for details. 

discuss whether LBVs may explode in Sect. 7. 

2. WOLF-RAYET MASS-LOSS RATES AS A 
FUNCTION OF METAL CONTENT 

In recent years, it has become clear that long- 
duration gamma-ray bursts (GRBs) are associated 
with the explosion of a mass ive star, providing impe 



tusto the collapsar model (jMacFadven fc Wooslev 
1999). The model works best if the progenitor fulfils 



the following two criteria: (i) the absence of a thick 
hydrogen envelope (enabling the jet to emerge) , and 
(ii) rapid rotation of the core (allowing a disk to 
form). This may point towards a rapidly rotating 
WR star. 

WR stars are believed to be the result of mass- 
loss durin g earlie r evolu tionary phases (the "Conti" 



scenario 
tary scenario, 



(IConti I Il976l) ). while in a complemen- 
the removal of the thick hydro- 
gen envelope may be due to a companion. Re- 
cently, an alternative scenario for p roducing a GRB 



progenitor has gained popula rity (|Yoon fc Langer 



120051 : 1 Wooslev fc Heger Il2006h : when a star rotates 
rapidly, it may mix "quasi homogeneously" , and the 
object may not develop the classical core-envelope 
structure, but remain small. A potential problem 
for producing a GRB in this scenario is that Galac- 
tic WR stars have strong stell ar winds whic h may 
remove the angular momentum ( Langer 1|l998? ), mak- 
ing it challenging to produce a GRB. 

This problem might be overcome if WR winds are 
weaker at low Z, so the question is: "are the winds 
of WR stars Z-dependent?" and if so, "how strong is 
this dependence?" The dense winds of WR stars are 
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Fig. 2. Wind efficiency 77 = (Mv 00 )/(L»/c) as a function 
of effective temper ature. The predictions are taken from 
IVink et al. I (|2000l ). Note the presence of the bi-stability 
jump around 25 kK, where 77 increases as Fe recombincs 
to Fe in. 



likely driven by radiation pressur e (jNugis fe Lamers 
2002uGrafener fe Hamann 1120051 ). just like their less 
extreme O star counterparts. 

This need not imply that WR winds depend on 
metal content, as WR stars produce copious amounts 
of metals such as carbon (in WC stars). If, on the 
one hand, these self-enriched elements dominate the 
driving (by their sheer number of particles), one 
would expect WR winds to be independent of their 
initial Z and the requirements of the collapsar model 
may never be met. If, on the other hand, iron (Fe) 
is pre dominantly respon sible for the driving (as in O 



stars; IVink et al. 1120011 ), WR winds might indeed be 
less efficient in low Z galaxies. 

To address t his issue regarding the .^- dependence 
of WR winds, (jVink fe de Koter 1 l2005f) computed 
mass-loss rates for late-type WN and WC stars as 
a function of the initial metal content (representa- 
tive of the host galaxy Z). The results are shown in 
Fig. [TJ For a discussion of the flattening in the mass- 
loss- dependence for initial metallicities below log 
(Z/ Zq) = — 2 and potential consequen ces for the firs t 
stars (Pop in), the reader is referred to I Vink ( 20061 ). 
but for the Z range down to log (Z/Zq) = —2, the 
mass loss is found to drop steeply, as M oc Z 85 , for 
the WN phase - where WR stars spend most of their 
time. This inefficiency of WR mass loss at subsolar 
Z may prevent the loss of stellar angular momentum, 
and may provide a boost to the collapsar model. 

3. MASS LOSS FROM OB STARS: ABSOLUTE 
RATES AND THE BI-STABILITY JUMP 

We now switch from a discussion of Z-dependent 
mass loss to one of T e ff-dependent mass loss. We de- 



scribe the expected wind properties in terms of their 
wind efficiency number rj = (Mvoo) I (L* / c), a mea- 
sure for the momentu m transfer from the photons to 
the ions in the wind. IVink et al. I (|2000h computed 
wind models as a function of effective temperature 
(Fig. [2]) . The overall behaviour is one of decreasing 
77 with decreasing T e g due to a growing mismatch 
between the wavelengths of the maximum opacity 
(in the UV) and the flux (gradually moving towards 
longer wavelengths). The behaviour changes at the 
"bi-stability jump" (BSJ; e.g. Lamers et al. 1995), 
where 77 increa ses by a factor of 2 -3, as Fe iv recom- 



bines to Fe ill (jVink et al. Ill999l b 



R ecent mass-loss studie s (jTrundle fe Lennon" 



20051 ; ICrowther et al. 1 120061 ) have reconfirmed dis- 
crepancies between empirical mass-loss rates an d 
predictions for B supergiants (jVink et al. ] l200fj . 



Discrepancies have also been reported for O stars 
(jBouret et al. 1 120031; iFullerton et al~ll2006l ). and it 
is as yet unclear whether the reported discrepancies 
for B supergiants are due to model assumptions (e.g. 
the neglect of wind clumping) or the physical reality 
of the BSJ. The most accurate way to derive M is 
belie ved to be through radi o observations. Intrigu- 
inglv. iBenaglia et al. ( 2007 ) present empirical radio 
mass-loss rates as a function of effective temperature 
that rese mble the mas s doss e fficiency behaviour pre- 
dicted bv lVink et all (|2000h . This may well be the 
first evidence of the presence of a mass-loss BSJ at 
the boundary between O and B supergiants. The 
relevance for stellar evolution is that when massive 
stars evolve at constant luminosity towards lower 
T e g, they are anticipated to cross the BSJ. Inter- 
estingly, LBVs brighter than log (L/Lq) = 5.8 (see 
Fig. [3]). are expected to encounter it continuously - 
on timescales of their photometric S Doradus vari- 
ability, discussed in the next section. 

4. MASS LOSS FROM LUMINOUS BLUE 
VARIABLES 

LBVs are unstable massive stars in the up- 



per p art of the HRD (e.g. iHumphrevs fe Davidson 
19941 ). As can be seen in Fig. [3l the "classical" 
LBVs, like AG Car, are anticipated to cross the BSJ 
at ~ 21 000 K. One of the defining characteristics 
for LBVs is their S Doradus (SD) variation of ~1 
- 2 mag on timescales of ye ars (short SD phases ) 
to decades (long SD pha ses) ( van Genderen Il2~00l[ ). 
Vink fe de Koter I (j2002h computed LBV mass-loss 
rates as a function of T e R- shown in Fig.JH Overplot- 
ted are the empiric al Ha mass-loss rates for AG Car 



(jStahl et al. 1120 01). which vary on the timescales 



of the photometric S Doradus variability. Although 
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Fig. 3. The LBVs in the HRD. The shaded areas rep- 
resent the S Doradus instability strip (diagonal) and the 
position of the LBVs during outburst (vertical). The 
dashed vertical line at 21 000 K indicates the posi- 
tion of the b i -stabi lity jump. The figure is taken from 
ISmith et al.l (|2004T l. 



the ag reement is not perfect (see IVink fe de Koter 
( 20021 ) for a discussion), the amplitude of the pre- 



dicted variability fits the observations well, and most 
importantly the overall behaviour appears to be very 
similar, and may indeed be explained in terms of 
the physics of the BSJ. T his bi-stable behaviour in 
an individual stellar wind ( Pauldrach fe Puis 199Clf ) 



causes the star to flip back and forth between two 
states: that of a low mass loss, high-velocity wind, 
to a high mass-loss, low velocity wind. The wind 
density (oc M/voo) would therefore be expected to 
change by a factor of ~2 x ~2, i.e. ~4 on the 
timescale of the SD variations. In the absence of any 
other material around the star, this would result in 
a pattern of concentric shells of varying density. 

5. RADIO SUPERNOVAE AND PROGENITOR 
MASS LOSS 

Radio SNe (RSNe) lightcurves and the model for 
SN interaction with the surr ounding circumstella r 
material has been reviewed by IWeiler et al.l (Il986l ). 
The radio emission is due to non-thermal electrons, 
while the absorption may be due to both syn- 
chrot ron self absorption as well as free-free absorr . 
tion ( Chevalier 19821 : Fransson fc Biornsson 1998 



Examples of the rise, peak, and power-law decline of 
radio lightcurves are shown in Fig. [5] (The episodic 
bumps at late time are discussed in Sect. [6]) 

The model constrains the wind density and thus 
the ratio of M to the termina l wind velocity (vm) : 
p oc M/voor 2 . Assuming Voo, IWeiler et al. I (2002) 
list M values in the range 10~ 6 -10~ 4 M Q yr _1 . For- 
tunately, these values agree with mass-loss predic- 
tions, but are broadly representative for massive 
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Fig. 4. Predicted (dotted line) and empirical (dashed 
line) mass-loss rates versus T e ft for the LBV AG Car. 
Note that both the qualitative behaviour and the ampli- 
tude of the mass-loss var iations are well reproduced. See 
IVink fc de Koter I i|2002f ) for details. 

stars over almost all post-main sequence evolution- 
ary phases, making it hard to infer the progenitor 
from radio lightcurves alone, unless these lightcurves 
betray their progenitor in some another way. 

6. QUASI-PERIODIC OSCILLATIONS IN 

RADIO SNE LIGHTCURVES 
A number of recent RSNe have shown sinusoidal 
modulations in their radio lightcu rves, in particu- 



lar SN 2001ig dRvder et al. 1 120041) and SN 2003bg 



(ISoderberg et al. I l2006f) ~ are strikingly similar in 



terms of both amplitude and variability timescale 
(see Fig. [5]). The recurrence timescale t of the bumps 
is ~ 150 days. Using Eq. (13) from (jWeiler et al 



198G): 



AP = 



R 



shell 



ejecta t'l 



id m 



(1) 



where m is the deceleration parameter (here m = 
0.85) and t[ is the time of measurement of the ejecta 
velocity relative to the moment of the explosion. As- 
suming u W ind = 10-20 km sec" 1 , typical wind veloci- 



ties for red (super)giants, (jRvder et al. 1120041 ) found 



a period P between successive mass-loss phases that 
was too long for red (super)giant pulsat ions (100s 



of days, see however iHeger et al 



( 19971 )). but too 
short for thermal pulses (10 2 -10 3 years). They 
therefore invoked an edge-on, eccentric binary sce- 
nario involving a WR-star and a massive compan- 
ion. One of the main differences between LBV and 
red giant winds is that LBV winds are about 10 
times faster. If the progenitor of SN 2001ig were an 
LBV, the expected period between successive mass- 
loss episodes would be AP ~25yr (for an assumed 
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Fig. 5. Radio luminosity versus time for two strik- 
ingly similar recent SNe: 2001ig and 2003bg. Note the 
quasi-sinusoidal modula tions during the power- law de- 
cline phase. Taken from ISoderberg et al. I (|2006h . 

200kmsec~ 1 ), cons istent with the long SD 



ph ase (IKotak fc Vink 

Soderberg et al. ( 20061 ) infer density enhance- 



ments of a factor of ~ 2 during the deviations from 
pure power-law evolution. They consider a range of 
options that might account for the modulations, but 
they favour a single-star progenitor model of a WR 
star that underwent episodes of intensified mass loss. 
However, they do not specify the physical mechanism 
that gives rise to these periods of enhanced mass 
loss. Our SD mechanism for LBVs may alleviate 
this shortcoming. 

7. DISCUSSION: DO LBVS EXPLODE? 

Are LBVs viable SNe progenitors? It may be rel- 
evant that both SNe 2001ig and 2003bg are "transi- 
tional" objects. SN 2001ig was initially classified as 
type II (showing H lines) but metamorphosed into a 
type Ib/c object (no H lines, weak He lines) about 
9 months later. This suggests that it has lost most 
of its H-rich envelope. SN 2003bg however was first 
classified as a type Ic, but within a month the spec- 
trum evolved into a type II SN. This transitional be- 
haviour hints at the fact that their progenitors are 
intermediate evolutionary objects: H-rich compared 
to OB/red (supcr)giants, but H-poor compared to 
WR stars. LBVs are likely candidates. 

Recently there has been much discussion regard- 
ing clumping in the winds of O stars. The value 
for the clumping fact or is very much an open issue. 



Mokiem et al. I (|2007h show that if the empirical Ha 
rates are overestimated by a factor of two due to 
clumping, these empirical rates are in good agree- 



(|2000l . 1200 ll ). and consequently our current knowl- 
edge of massive star evolution is not anticipated to 
be affected by clumped winds. If however the wind 
clumping factor would be significantly larger than a 
factor two/three (as has been suggested by UV anal- 
yses) , this could have severe implications for massive 
star evolution. One consequence might be that giant 
LBV eruptions (j] Car type eruptions, not the typify- 
ing SD variations) dominate the integrated mass loss 
during evolution ( Smith fc Owocki l200ot ). An al- 



ternative scenario could be that post-main sequence 
stars do not become WR stars, but explode early - 
during their LBV phase. 

Here, we have presented indications that at least 
those SNe that show quasi-periodic modulations in 
their radio lightcurves might have LBV progenitors 



(|Kotak fc Vinkll2006[) . It has also been speculated 



that LB Vs may be the generic progenitors of type 
Iln SNe ( Gal- Yam et al. 2006h . however it may be 
more relevant to discuss Iln SNe as a "phenomenon" 
describing SN ejecta expanding into a dense CSM 
than a one-to-one correla tion to a particular progen- 
itor (jKotak et al. 1 120041 ). Nevertheless, some frac- 



ment with the mass- loss predictions of I Vink et al 



tion of type Iln SNe may well have LBV progenitors 
although the observational evidence remains elusive. 

It is relevant to note that the LBV candidate 
HD168625 is embedded in a bipolar-shaped neb- 
ula that resembles the triple-ring system around 
SN1987A. This similarity could hint that the pro- 
genitor of 1987A (i.e. the blue supcrgiant Sk-69 202) 
underwent an LBV giant eruption before it exploded 
(ISmithll2007h . 

Future mass-loss predictions are anticipated to 
play an important role in obtaining knowledge about 
the lives and deaths of massive stars. 
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